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Abstract 

We analyze the Charlier polynomials C„(x) and their zeros asymptotically as n ^ 
oo. We obtain asymptotic approximations, using the limit relation between the Krawtchouk 
and Charlier polynomials, involving some special functions. We give numerical exam- 
ples showing the accuracy of our formulas. 
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1 Introduction 



The Charlier polynomials C„(x) ^ are defined by 



CJx] 



2-^0 



-ra, —X 



where x>0, n = 0,l,... and a > 0. They satisfy the discrete orthogonality condition |37j 
They are part of the Askey-scheme ^H] of hypergeometric orthogonal polynomials: 
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i \ 
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/ i \ / i 
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i 
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Laguerre 



Charlier 



2-^0 



Hermite 



where the arrows indicate limit relations between the polynomials. 

The Charlier polynomials have applications in quantum mechanics j2Hl, IHIji P^ - 
difference equations [231 , teletraffic theory [TB], j27j, generating functions [22], [20], 
and probability theory [3], [21], (HO], [^- The g-analogue of the Charlier polynomials were 
studied in [2], |H], JH] and [H]. The generalized Charlier polynomials were analyzed in |14j . 
dl], jSl, jSl] and [10]. 

Asymptotics for the L^-norms and information entropies of Charlier polynomials were 
derived in [21]. Bounds for their zeros were obtained in [20]- Asymptotic representations 
were established in ^1] in terms of Hermite polynomials and in [21] in terms of Gamma 
functions. Some asymptotic estimates were computed in from a representation of C„(x) 
in terms of Bell polynomials. An asymptotic formula when x < was derived in [2S] using 
probabilistic methods. 

In [13], Goh studied the asymptotic behavior of C„(a;) for large n using an approximation 
of the Plancharel-Rotach type. A uniform asymptotic expansion was derived in [B] using the 
saddle-point method. Asymptotic expansions were obtained in [TIT] from a second order linear 
differential satisfied by in which a is the independent variable and a; is a parameter. 
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In this paper we shall take a different approach and investigate the asymptotic behavior 
of Cn{x) as n oo, by using the limit relation between the Krawtchouk polynomials Kn{x) 
defined by 



—n, —X 
-N 



Kn{x)=Kn{x,p,N) = 2F1 

and the Charlier polynomials, namely 



n = 0, 1, . . . , A^, < X < iV, < j9 < 1 (2) 



limKr^ix,^ N)=Cn{x). (3) 

iV— >oo V iV / 

We shall use the asymptotic expansions derived in ^ for the scaled Krawtchouk polynomials 
kn{x), with 

kn{x) = kn{x,p,N) = {-pY (^^^Kn{x,p,N). (4) 

A similar idea has been used in J2] and [38j to obtain asymptotic approximations of 
several orthogonal polynomials of the Askey-scheme in terms of Hermite and Laguerre poly- 
nomials. 



2 Preliminaries 

The following is the main result derived in jH]. 

Theorem 1 As N 00, kn{x,p, N) admits the following asymptotic approximations (see 
Figur^. 



1. n = 0(1), < ?/ < 1, y 



where 



kn{x)^k^^\y) = —{y-pr (5) 



y z 
e = N^-^, X = -, n = - < y, z < 1. 

e e 



2. n = 0(1), y ^ p, y = p + riy/2pqe, rj = 0(1). 

fcn(x)~A:f(r/) = ^(f)^H.(r^), (6) 
where q = 1 — p and H„ [rj) is the Hermite polynomial. 
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Figure 1: A sketch of the different asymptotic regions for kn{x). 
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3. <y < Y~{z), < z <p, where 



with 



Y^{z)^p+{q-p)z± 2zUo, Uo{z) 



pq{l - z) 



kn{x) - k^'\y,z) = ^exp [i;{y, z,U-)e-'] L{z,U-), 



z)^{z- 1) ln(C/) + (1 - y) HU -p)+y ln(C/ + q), 



L{y,z) 



l{ U-p){U + q) 
[U' - {Uof] 



and 



2 \ z 

4. Y+{z) <y<l, 0<z<q. 

Ux) ~ k^'\y,z) = ^exp [il;{y,z,U+)e-'] L{z,U+). 



27r 



5. X ^ 0(1), p < z <1. 



kn{x) ~ k^^^ {x, z) — , ^ cosifix^ ^ ^ 



2'K^Z{\ - z) 



V 



exp [0o(^)£ ^] 



T(x) sin(7rx) ( — — j exp 



[z — 1) ln(g) + T^\z 



TT Z — p \Z — p ^ 

where 

4>o{z) = — 1) ln(l — z) — ^ln(^) + zln{—p) 
and r{x) is the Gamma function. 

6. X — z p, z — p — Uy/pqe, u — 0{1). 



*„M~ 




X exp 



Trip — gin (g) Uy/pqKi—Uy/pqhi^q) u 



21 



where Dx{u) is the parabolic cylinder function. 
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COS 



7. < y < Y-{z), p< z <1. 

kn{x) - k^'^\y, z) = exp ^ 

8. Y-{z), 0<z<p,y = Y-{z) - (3e^/\ /3 = 0(1). 

kn{x) ~ k^^\P, z) = exp 
where 



(^"^ k('\y,z) + 2isin (^) k'^'\y,z)\ (15) 



Ai(e'p)-^, (16) 



M^) = -^TTi + (^ - 1) In (Uo) + Y-{z) In (C/q - q) + [l - Y'iz)] In (C/q + p) , (17) 



anc? Ai (•) is the Airy function. 




z (Uo+p)(Uo-q) 



(18) 



P. y Y-{z), p<z<l, y = Y-{z) (3 = 0(1). 

^„(x) ~ A;(^)(/3, z) = £^ exp L(^)£-' + In (^f^) /^^"^ 

X [A+(/3,z)Ai +iA-(/3,^)Bi (^^/j)] , 

where 't}{z) = —Q{z), 

A^(/3, ^) = exp |27ri (z) - /^si s"^} ± 1. 
anc? Ai (■) , Bi (■) are the Airy functions. 
10. Y-{z) <y< Y+{z), 0<z<l. 

kn{x) ~ k^''\y.z) = k^''\y,z)-rk^^\y,z). 



(19) 



(20) 



(21) 



11. y K,Y^{z), Q < z < q, y = Y+{z) + ae^, a = 0(1). 

= £3 exp 



V'i(^)£-' + ln(^^la£ 5 



Ai 



(ei)3 a 



(01 



(22) 



where 



^P^{z) = {z-l) In (t/o) + Y+{z) In (^/q + g) + [l - Y+{z)\ In (^7o - p) , (23) 



ei(^) = 



z {Uo-p){Uo + q) 



(24) 



In order to obtain the corresponding asymptotic expansions for Kn{x), we need to derive 
asymptotic formulas for (— p)" (^) in the different regions of Theorem 1. 
The following lemma follows immediately from Stirling's formula jl] 

r(x) ~ W— x"e-", x^oo. (25) 
V X 

Lemma 2 As N oo, we have the following asymptotic approximations: 



'-pT i^) - \ n = 0(1). (26) 



'-VT - f' , exp [<P{z)e-'] , n = ze-' (27) 

y l-n^ z (1 — z) 



n 



where 



3. 



where 



(f){z) = z ln(p) + zTTi + {z-l) ln(l - z) - z \n{z). (28) 
;-j9)"f^') ~-^exp[0i(«)], n = pe-'-uM, u = 0(1) (29) 
4>i{u) = [nip — q ln(g)] e^^ — u^ppq [ln(g) + 7ri] £"2 — -v?. (30) 



3 Limit analysis 

Setting 

y = xe, z = ne, q = 1 — p, p = ae (31) 
in (|7j) and letting e — > 0, we obtain 

Y^{z)^X^{a,n) = {V^±V^)\ U,{z) ^ (32) 

Hence, the eleven regions of Theorem 1 transform into the following regions (see Figure Ej)- 

1. Region I 

From (jni) and we have for n = 0(1) 

Kn{x) ~ (1 - ^) • (33) 
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Figure 2: A sketch of the different asymptotic regions for C„(x). 
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Thus, 



C„(x) 



a, 



The formula above is exact for n = 0, 1 and in the hmit as a — > cxd we have 



hm Cn{av) = (1 - vY = iFq 



-n 



2. Region II 

From © and (j^Bj) we have for x = pe~^ + r]\ rj = 0(1) 



^n(x) ~ (-1)" ( I ) H„(r/). 



Hence, 



~ (-l)"(2a)"2 H„(ry), x = a + r/V2a. 
Equation ()36|) is exact for n = 0,1 and in the limit as a ^ oo we have 

lim (-1)" (2a) 2 C„(a + r/v^) = H„ (r/) 

a— >oo 

which is equation 2.12.1 in. 
3. Region III 

From dHl), (EZI) and ^ we have for < ?/ < Y-{z), <z <p 



Kn{x) ~ K^^\y,z) = exp 

where 

0(2, t/) = 

From (0), (dH) and dSHl) we obtain 



{l-z){U-p){U + q) 

[/2 - [/o2 



G(2,f/-) 



'^iy,z,u ) -0(2) 



^3(2;) 



where 



, . X , ,'a + x-ra + A\ , fa-x + n + A\ 1. 

W^lxl = X m + n in H — [a — x — n 

^ * 2a / V 2a 2^ 



and 



G(z,f/-) ^L3(x) = 
10 



a — X — n + A 



2A 



X 



Figure 3: A comparison of C„(x) (solid curve) and -^3(0;) (000) for n = 30 with a = 50.165184. 

for < X < X~ , < n < a, with 

A{a,n,x) = \/ — 2a{x + n) + {x — n)^. (41) 

Thus, 

Cn{x) ~ Fsix) = exp [^3(x)] Ls{x), 0<x< X', < n< a. (42) 

Figure El shows the accuracy of the approximation (j42|l with n = 30 and a = 50.165184 
in the range —3 < x < X~ . 

Remark 3 Although we have not shown proof of it, the approximation is in fact 
also valid for x < and n > (see Figure^. 

4. Region IV 
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Figure 4: A comparison of C„(x) (solid curve) and F^{x) (ooo) for n — 30 with a — 2.165184. 
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From (O, (EZ|) and ^ we have for Y+{z) < y < I, < z <q 



Kn{x) ~ K^'^\y, z) = exp 



'ij{y,z,U+)-4>{z) 



G(z,f/+). 



(43) 



From ©, (HH) and dnni) we obtain 

ilj{y,z,U+)-<P{z) 



\&4(x) — nvri, 



where 



^ , , , ,'a + x — n — A\ , /x — a — r2 + A\ 1, 

^4(x) = X In ( + n In I +-{a-x-n + A) (44) 



and 



G{z,U-) ^^4(5 



X — a + 72 + A 



2A 



for < X. Therefore, 

Cn{x) ~ F4(x) = (-l)"exp [^4{x)] U{x), X+ < X. 



(45) 



(46) 



Figure El shows the accuracy of the approximation ()46p with n = 30 and a = 2.165184 
in the range < x < 00. 

5. Region V 

From (0, (EZ|) and ((211) we have for x = 0(1), p< z <1 



Kn{x) ~ exp 



xln ( 1 

P 



Z — p \ TC 



^ / 2 

cos (TTx) ■ — \l —z (1 — 2;)r (x + 1) sin (ttx) (47) 



X exp 



(l-^)ln(i^)+zln, , , 

+ X m 



z — p 



Hence, 



Cnix) ~ F5(x) = exp 



X In ( 1 

a 



cos (ttx) — v^a/ — r (a^ + 1) sin (ttx) (48) 



X exp 



n In 



— ) — (x + 1) ln(n — a) + a — n 
a/ 



X ~ 0, n > a. 



Figure ini shows the accuracy of the approximation PH|) with n = 30, a = 2.165184 and 
X ^ 0. 
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Figure 5: A comparison of C„(x) (solid curve) and F4{x) (ooo) for n — 30 with a — 2.165184. 
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Figure 6: A comparison of C„(x) (solid curve) and ^5(2;) (000) for n — 30 with a — 2.165184. 
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Figure 7: A comparison of C„(x) (solid curve) and Fq{x) (ooo) for n = 30 with a = 30.165184. 
6. Region VI 

From (fT^ and (^^ - (jHnj) we have for x = 0(1), z = p — u^/pqe, u = 0(1) 

,2" 



Knix) ~ exp 



11] + _ 



(u) 



(49) 



Thus, 



Cn{x) ~ F6(x) = exp 



Da; (m) , x ~ 0, n = a — u\fa. (50) 



Figure [7| shows the accuracy of the approximation (jSUj) with n = 30, a = 30.165184 
in ~ a) and x ~ 0. 



7. Region VII 
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From (Uni), (|SZ|) and (gSI) we have for < y < Y-{z), p< z <1 



K„ix) ~ exp l^^j [cos (^) z) + 2isin (^) ^^(^^(y, ^) 

Therefore 

C'n(a^) ~ exp (vrix) [cos (ttx) C^\x) + 2i sin (7ra;)C(^)(x)] 
for ^ X < X~ , n > a, which we can rewrite as 



(51) 



Cn{x) ~ Fj^x) = exp 



, n — a — X + A\ / a + n — X — A\ 1 

X m + nm \ H — (a — x — n 

2a / V 2a / 2 ^ 



/x + ra — a + A , /x + n — a — A 

X cos iiTx] \l : 2 sin ( nx 



X exp 



2A 

, , n — a — X — A . . 

X In I J + n m 

2a 



2A 

a + n — X + A\ 1 



(52) 



2a 



- (a — X — n + A) 



Figure IHl shows the accuracy of the approximation (j52j) with n = 30 and a = 2.165184 
in the range <^ x < X~. 

Region VIII 

From (Uni), (EZ|) and (j2Hl) we have for y ^ Y-{z), < z < p, y = Y-{z) - /Je^/s^ 
/3 = 0(1) 



Kn{x) ~ £ 6 exp 



e ^\Uo-qJ 



(53) 



X \/27r 



z{l-z) 



pq 



Ai e'/? e-3. 



Since 

we have from ()32j] 
From (Hlj), (I2HI) and dnni) we get 



P= [Y-{z)-y] s- 
I3e-^ ^X- -X. 



(54) 
(55) 



m^im + In ( ^) pe-'^ ^ 1. In (!^) + xln ( 1 - . 1 + - v/^. 
e \Uq — qj 2 \a/ V Va 

From ()18|1 and we obtain 



e 6^277 



(56) 
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Figure 8: A sketch of Cn{x) / F-j{x) for n = 30 with a — 2.165184. The vertical hnes are due 
to the discontinuities at a; ~ 13, 14 and 15. 
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and 



Therefore, 



n\s (X — x) 



2 • 

3 



(57) 



Cn{x) ~ F8(x) = ( -) ' (v^ - v^) 5 Ai 



n\i (X — x) 



(58) 



X exp 



-nln f — ^ + x\n ( 1 
2 Va/ V 



n 
a 



for x~X , 0<ri<a. 
9. Region IX 

From (Uni), (1201), (I2Z|) and ^ we have for y ^ Y-{z), p < z < 1, y = Y-{z) - 
P = 0{1) 



Kn{x) ~ e e exp 



il)o{z) - (t){z) fUo+p\ 1 



+ In 



X V27r 



f/o-g 

A+(/5, ;z)Ai + iA-(/3, z)Bi (^t?*/? 



which can be written as 



Kn{x) e 6V27rexp 



\z - 1) In {^) + F- In (g - Uq) + (1 - F") In {p + Uq) + z In 



(59) 



X exp 



In 



q-UoJ 



z{l-z) 



pq 



^-3 



cos(7rx)Ai h9'/5 -sin(7rx)Bi h?''/? 



Using and (|37|) in with = —9, we have 



Cn{x) ~ -F9(x) = v^27r ( — (i/n — i/a) ^ exp 

X "I cos (ttx) Ai 
for X ~ X~, n > a. 



n 



-nln— ) + xln \ 1 + Jan — n 



(60) 



n\l (X — x) 



sin(7rx)Bi 



n\i (X — x) 
«^ (v^-v^)^J 



19 



-2e+24 



-4e+24 - 



-6e+24 



-8e+24 - 



-1e+25 



-1 .2e+25 



31 




Figure 9: A comparison of C„(x) (solid curve) and -Fio(a;) (ooo) for n = 30 with a = 2.165184. 



10. Region X 

From dn}, dSZI) and (0 we have for Y'^z) <y < Y+{z), < ^ < 1 



Thus, 



(61) 



(62) 



Figure ini shows the accuracy of the approximation with n = 30, a = 2.165184 in 
the range X~ < x < X^. 

11. Region XI 



From (1121), (EZ|) and (UHl) we have ior y ^ Y+{z), < z < q, y = Y+{z) + as 



3 , a 



20 



0(1), 



Kn{x) ~ e 6 exp 



ipiiz) -(j){z) fUo + q \ _i 
h m I — I ae ^ 



X V27r 



^(1-^) 
pq 



Ai 



Since 

we have from 
From 



a= [y-Y+iz)] e 



ae 'A X — 



and we get 



e \Uo — p J 2 Va/ \ yet 



an — 



From ()24|) and ()65|) we obtain 



^(1-^) 
pq 



e 6 V27r 



(01 



/2^(^)'(v^+v^) 



and 



Therefore, 



(61)3 a 



2 • 

3 



Cn{x) ~ Fn(x) = (^^) ' (v^ + v^) ^ Ai 
X (— l)"exp 



n\-e (x-X+) 
«^ (v^+v^)^ 



-72 In f — ) + xln ( 1 

2 \a 



an — \/n 



for a; ~ X 



Comparison with previous results 

'^e shall now compare our results with those obtained previously in ^ and 

1. Region VII: < x < X~, n> a. 
Setting X = un, with 

^ / X la a 

u = 0(l), 0<M< 1-2W- + - < 1 

V n n 



21 



in we have, as n oo 



cosiunn) 
Frix) ~ griu) = exp 



\ uln(-)+{u-l)\n{l-u)-u\n + -^\ (69) 



-2 sin (unir) 



u 



l-u 



exp 



In 



n 



+ (1 - m) In (1 - m) + u ln(u) - 1 



n + 



l-u 



The second term of equation (1^^) is the same as the equation before (5.3) in and 
equation (84) in ^^l- However, the first term is absent in previous works, although it 
is necessary in the asymptotic approximation, especially when m ~ 0, 1, 2, . . . . 

2. Region IX: x ~ X~, n > a. 

We now set x = X~ + tn^ , t = 0(1) in (jUUI) and obtain, as n — >■ cxo 



Fg{x) ~ gg{t) = V2na exp 



^ + tn^ + In ^ 



X < cos 



X- + tnTi 



TT 



Ai f —ta 6 ) — sin 



X" + trie ] IT 



3 

n H — a 
2 

Bi f -ta~t 



Equation (fTUI) agrees with equation (5.13) in P| and equation (51) in [T3] . 

3. Region X: < x < X+. 

Setting a; = n + a + 2 sin(6') v^on, with — | < < | in we have, as n — > cxd 



_ 1 1 

i^3(x) ~ 93{e) = (-1)" ^^exp 



ln(!^ 

a 



1 



TT 

n + -i 
4 



X 

X exp <^a 
Similarly from ()46p we get 

X exp < i/on 



v/2 cos (^) 

exp I ^/an sin (0) In - sin (^) (20 - tt) i - 2 cos (61) i > 
1 - i cos (20) + i In (^) - 1 sin (26) i-Oi + ^i 



_ 1 1 

a 



exp 



ln(!^ 

a 



IT 

n — -i 

4 



X 



^2 cos (0) 

sin (0) In (^^) + sin {&) (20 - tt) i + 2 cos (0) i } 
exp I a 1 - - cos (20) + - In [-] + - sin (20) i + 0i - -i I 



(70) 



(71) 



(72) 
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Using (HH) and ^ in ()62|) we have 



Fioix) ~ gioi0) = (-1) 



exp 



In 



n 



n 



X exp < V'^'T- 



sin {6) In 



a/ J 



+ a 



l-icos(2«) + ll„Q]} 



X cos I A/an [sin (6*) (26* - tt) + 2 cos (9)] + a 
Equation (|73|l is equivalent to equation (44) in [T3] . 



lsin(2^^) + e-^ 



4. Region XI: x X+. 

We now set x = + sn^, s = 0(1) in ()68p and obtain, as ^ oo 



- ( + sne + n] In f — ^ — n + -a 
2 V J \aJ 2 



-P'ii(a;) ~ 9n{s) = V2na en^ exp 
Equation (f7I|l is equation (5.12) in jB] and equation (30) in [T 



(73) 



Ai f sa 6 



(74) 



5 Zeros 

Using the formulas from the previous sections we can obtain approximations to the zeros of 
the Charlier polynomials. 

1. x ~ 0, n > a. 

The first zero is exponentially small. From ()48p we have, as x — > 
C^{x) ~ 1 



In ( - n - ^a-n-e'^- 



a 



n — a 



X. 



Solving for x we obtain 



Xo 



V2 



ln(^-l 

a / n — a 



V2 



n ^ oo 



Tin 



where xo denotes the smallest zero. 

2. < X < X~, n> a. 

In this range of x, the zeros are exponentially close to 1,2,..., [X~ \ . Using 

X — n — a + 



t 



1 



(75) 



23 



and the asymptotic formulas ^3] 



Ai(x) 



Bi(x) 



exp 



2 ^ 
-1X2 



2y/TTX4. 





2 i^' 


exp 


3X2^ 



X ^ 00 



a; ^ 00 



we have, as n — >■ 00 



Ai - 



-tae 



Bi - 



1 

-tae 



exp 



4n 4 (^X — x) 



Using (fTBj) in (fTOj) we have 



1 



g^it) ~ ^ - exp 



— a in 4 iX — x) 
3 ^ ' 



tan (ttx) . 



Since x^ ~ j, j = 1, 2, . . . , [X J , we get 

1 



2 exp 



-^a 4n 4 (X - j) 



7r(x. - jj 



which we can solve to obtain 



Xj ~ J H — exp 



3. X- < X < X+. 



(76) 



j = l,2,...,[X-J. (77) 



Finally, the non-trivial zeros of the Charlier polynomials can be approximated using 
dZHl). We have gio{e) = if and only if 



cos I Van [sin (6) {26 - vr) + 2 cos (9)] + a 
or equivalently if 

[sin (6) {26 - vr) + 2 cos (6)] + a 



- sin (2^) + 6-- 







1 sin (2^^) + - I 



vr TT , , ^ 
- = - + vr/, I ez 
4 2 



or 



/an [sin (6) {26 - vr) + 2 cos (6)] + a 



lsin(2^) + 0-^ 



37r 



nl = 0, (78) 
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with -| <6 < |. Recalling that 



X = n + a + 2 sin{0)^/(m, 



(79) 



we see that the condition X < x < X+ implies 



< / < 2J(m-a - -. 

_ _ V ^ 



(80) 



Equation (fTSj) cannot be solved exactly. However, it can be easily solved numerically to 
any desired accuracy and using (fT^ll gives very good approximations for the nontrivial 
zeros. 

In Table 1 we computed the exact and approximate zeros of 6*25(0;) with a = 2.16564899 
using (HSl), (HID and ((ISll-dHni)- 

Conclusion 4 We analyzed the asymptotic behavior of the Charlier polynomials in the range 
< X as n ^ 00. We also obtained approximations for their zeros. We intend to extend 
our method to the other polynomials of the Askey-scheme to obtain asymptotic expansions of 
them. 
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